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E(A+M)PEC - An OpenCL Atomic & Molecular Plasma Emission 
Code For Interstellar Medium Simulations 
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Abstract. E(A+M)PEC traces the ionization structure, cooling and emission spectra 
of plasmas. It is written in OpenCL, runs in NVIDIA Graphics Processor Units and can 
be coupled to any HD or MHD code to follow the dynamical and thermal evolution of 
any plasma in, e.g., the interstellar medium (ISM). 



1. Introduction 

The thermal evolution of the ISM is determined by heating and cooling processes which 
may not be synchronized with ionization and recombination processes and the system 
deviates from coUisional ionization equilibrium (CIE). A step forward in ISM mod- 
elling requires that both the thermal and dynamical evolutions be treated in a self- 
consistent way and under non-equilibrium ionization conditions (NEI). 
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Figure 1 . Normalized to n^ NEI (left: solid lines and dashed lines show cooling 
under isochoric and isobaric coditions, respectively) and CIE (right) cooling effi- 
ciencies using I Asplund et al. ( 2009) solar abundances. Total cooling (black lines), 
free-free (brown line), radiative (red) and dielectronic (green) recombination, elec- 
tron impact ionization (blue), two-photon {2y; orange) and line excitation excluding 
2y (magenta) emissions. 
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2. E(A+M)PEC Features and Some Results 

E(A+M)PEC traces the ten most abundant elements in nature, includes up to date solar 



abundances and atomic and molecular data (see e.g., http://www.lca.uevora.pt/res earchl 
and references therein). The physical processes comprise electron impact ionization, 
excitation auto-ionization, radiative and dielectronic recombination, charge-exchange 
reactions, ionization of Hi as result of He ions recombination, continuum (free-free, 
free-bound, two-photon) and line mission (A e [lA, 610ju]), and H, C and O chemistry. 
The cooling efficiencies of a plasma evolving from 10^ K (where it was fully 
ionized) to 100 K under CIE and NEl (isochoric and isobaric) conditions are shown in 
Figure [I] Departure from CIE occurs at some 10^ K. As recombination lags behind, 
single and double ionized species exist at lower temperatures (Figure |2]l something 
that does not occur under CIE (because ionization and recombination is synchronized) 
and therefore neutrals form at temperatures near 10"^ K. These differences also affect 
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Figure 2. CIE (solid lines) and NEI isochoric (dashed lines) ionization structure 
of Carbon between 10^ and 10^ K. Charge-exchange reactions have been included. 

the emission spectra - with recombination lagging behind, typically around 10^ K, the 
spectra at A < 100 A becomes free-bound dominated under NEl conditions. 



3. Final Remarks 



The ionization structure, cooling and emission spectra obtained with the E(A-i-M)PEC 
running in NVIDIA CPUs are similar (though with some deviations due to different 



atomic data and abundances) to results of, e.g, Schmutzler & Tschamuter (1993 1, Gnat 



|& Sternberg' (2007) and |Bryans e t al.'(2009). Regarding ISM simulations, it should be 
stressed that the interstellar cooling function is changing in space and time (Avillez & 
Breitschwerdt, submitted). 
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